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We report the observation of the decay B~ —> D^+ n - ^oK' , where D^+^-^o indicates a neutral 
D meson detected in the final state 7r + 7r _ 7r°, excluding Kgir . This doubly Cabibbo-suppressed 
decay chain can be used to measure the CKM phase 7. Using about 229 million e + e~ — * BB events 
recorded by the BABAR experiment at the PEP-II e + e~ storage ring, we measure the branching 
fraction B(B~ — > D w + n - w oK~) = (5.5±1.0 (stat.)±0.7 (syst.)) x 10 -6 and the decay rate asymmetry 
A = 0.02 ± 0.16 (stat.) ± 0.03 (syst.) for the full decay chain. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The Cabibbo-Kobayashi-Maskawa (CKM) matrix, 
whose element Vy describes the weak charged-current 
coupling between quark flavors i and j, provides an ex- 
planation for CP violation in the Standard Model. A 
crucial part of the program to study CP violation is the 
measurement of the angle 7 = arg {-V ud V* b /V cd V* b ) of 
the unitarity triangle related to the CKM matrix. The 
decays B -> D^°K^ can be used to measure 7 with 
essentially no hadronic uncertainties, making use of inter- 
ference between b — » ucs and b — > cus decay amplitudes. 
A number of variations on the original method 0] have 
been developed, and some have been explored experi- 
mentally. Employing multiple methods helps to resolve 
discrete ambiguities and decrease the experimental error. 

An important class of 7 measurement methods involves 
B -> with multi-body D decays 1 . In this 

technique, 7 is extracted from an analysis of the D- 
decay Dalitz plot, and ambiguities are resolved through 
interference between several D decay amplitudes [3j- 
Both Belle 4j and BABAR have used this method 
to obtain limits on 7 with the Cabibbo-favored decay 
D — > K®ir + ir~. The same approach can be carried out 
with multi-body final states that are produced by singly 
Cabibbo-suppressed decay of both D° and D° [|| . While 
these modes yield much smaller event samples, their in- 
terfering D° and D° decay amplitudes have similar mag- 
nitudes. Therefore, their overall sensitivity to 7 is a pri- 
ori expected to be similar to that of Cabibbo-favored D 
decays, where the interfering amplitudes typically have 
very dissimilar magnitudes. 

Among the singly Cabibbo-suppressed modes, the de- 
cay D — > 7r + 7T~7r° has a relatively large branching frac- 
tion 7] and a simple Dalitz plot dominated by broad 
p resonances making it attractive for the measure- 
ment of 7 with this technique. Its major difficulty is 
the relatively small signal-to-background ratio, which re- 
sults mainly from the high combinatorial background 



1 We use the symbol D to indicate any linear combination of a D 1 
and a D° meson state. 



associated with 7r° reconstruction. In this article, we 
describe an analysis procedure with which to extract 
the B~ — > D^+^-^oK" signal for later use in a Dalitz 
plot analysis measurement of 7, and report the mea- 
sured branching fraction and decay rate asymmetry of 
this decay chain. Our result excludes the decay mode 
D — > K®tt°, which is a previously-studied CP-eigenstate 
not related to the method of Ref. Q . 

The decay rate asymmetry A = (N + — N~ ) / (N + + 
N~), where N + (N~) is the number of signal B + (B~) 
decays, depends on the weak and strong phases of the B 
decay, as well as the D° and D° decay rate and phase 
variation over the Dalitz plot. Its magnitude is at most 
of order 2rs, where t\b, estimated to be about 0.1 0, 
is the ratio between the magnitudes of the interfering 
b — > ucs and b — > cus amplitudes. Due to interference, 
the branching fraction B(B~ — * D^+^-^oK - ) may differ 
from the product B pmd = B(B~ -> D°R-) x B{D° -> 
tt+tt-tt ) = (4.1±1.6)xlCT 6 7] by up to about 2r B B pmd . 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II energy-asymmetric e + e~ 
storage ring. The data consist of 207 fb _1 collected on 
the T(4S) resonance (on-resonance sample), and 21 fb _1 
collected at an e + e~ center-of-mass (CM) energy approx- 
imately 40 MeV below the resonance peak (off-resonance 
sample). Samples of simulated events were analyzed with 
the same reconstruction and analysis procedure. These 
include an e + e~ — > BB sample about three times larger 
than the data; a continuum e + e~ — > qq sample, where q 
represents & u, d, s, or c quark, with equivalent luminos- 
ity similar to that of the data; and a signal sample about 
200 times larger than what is expected in the data. The 
BABAR detector, as well as the methods used for charged 
and neutral particle reconstruction and identification are 
described in detail in Ref. 

We select events using criteria designed to maximize 
the signal branching fraction sensitivity and the reliabil- 
ity of the maximum likelihood fit procedure described be- 
low. To suppress the continuum background, we require 
the ratio H2/H0 of the 2nd to the 0th Fox- Wolfram mo- 
ments 10 . computed from the momenta of all charged 
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particles and photon candidates not matched to tracks, 
to be less than 0.50. Charged kaon candidates are re- 
quired to have a high quality particle identification mea- 
surement and be identified using kaon selection criteria 
that reduce the pion background to less than 3%. The 
measured energy of photon candidates is required to be 
at least 30 MeV. Photon candidate pairs whose invariant 
mass is within 25 MeV/ c 2 of the nominal 7r° mass [jj are 
combined to make ir° candidates, to which we perform 
a constrained-mass fit in order to improve the ir° energy 
and momentum resolutions. Throughout this article, we 
use the symbol to refer to the harder (higher-energy) 
of the two photons constituting a ir° candidate, and j s 
to denote the softer (lower-energy) photon. 

We select D — > tt + tt~tt° candidate decays by re- 
quiring the ir + ir~ir° invariant mass mo to be between 
1.830 GeV/c 2 and 1.895 GeV/c 2 . The m D resolution is 
about 14 MeV/c 2 . The D candidate energy and mo- 
mentum resolutions are then improved by performing 
a constrained mass fit. The charged pion candidates 
are required to fail kaon selection criteria. The decay 
D — > K®ir° is rejected by excluding iy + tt~ candidate 
pairs whose invariant mass is between 0.489 GeV/c 2 and 
0.508 GeV/c 2 . We note that this last requirement will 
not be needed when measuring 7 with an analysis of the 
7r + 7T~7r° Dalitz plot, where the K®ir° final state can be 
included as an incoherent term, as done in Ref. 

Candidate B~~ — > D n + n - n oK~ decays are constructed 
by combining a D — » 7r + 7r~7r candidate with a charged 
kaon candidate. Additional continuum suppression is ob- 
tained by requiring |cos#t| < 0.8, where 9t is the an- 
gle between the thrust axis calculated in the CM frame 
with the daughters of the B candidate and the thrust 
axis of the rest of the event (ROE). For each B can- 
didate we calculate the beam-energy substituted mass 
hues = ~ Ipb| 2 , where the total CM energy 

Eqm is continuously determined from the measured PEP- 
II beam energies, and pb is the momentum of the B can- 
didate in the CM frame. Signal events have a Gaussian 
rriES distribution that peaks at the nominal B~ mass 
with a width of about 2.7 MeV/c 2 , while background 
is distributed more broadly than signal. We require 
5.272 < m ES < 5.300 GeV/c 2 . The energy difference 
AE = E B - Scm/2, where E B is the CM energy of the 
B candidate, is required to be between —70 MeV and 
60 MeV. The AE distribution of signal events peaks 
around MeV with a width of 21 MeV. 

About 25% of the events selected have more than one 
B candidate. In these events, we select one B candidate 
at random. Random selection allows consistent studies 
of background suppression variables, and degrades the 
signal sensitivity by only a few percent relative to the 
best possible selection method. 

Studying the simulated event sample selected by the 
above criteria, we identify ten event types, one signal 
and nine background. We list these types with the labels 



used to refer to them throughout the article: 

• DKjj : B~ — > D^+^-^oK - events that were cor- 
rectly reconstructed. These are the only events con- 
sidered to be signal. 

• DKp : B~ — > D^+^-^oK - events in which the D 
candidate is misreconstructed, namely, some of the 
particles used to form the final state tt + tt~tt° do 
not originate in the decay of the D meson. 

• Dtt£> : B~ — > D°ir~ , D° — > ir + ir~ir° decays, where 
the decay D° — > 7r + 7r~7r° is correctly reconstructed 
and the remaining n~ is mistaken to be the kaon. 

• Dirp : B~ — > D°n~ , D° — > ir + ir~ir° decays, where 
the D candidate is misreconstructed. The kaon 
candidate may be either the remaining tt~ or a par- 
ticle from the other B meson in the event. 

• DKX : B -> excluding D -> tt+tt-tt 
decays, with a misreconstructed D candidate. 

• DttX : B -> DW%- and B -> D^p~, excluding 
D° — > 7r + 7T~7r° decays, with a misreconstructed D 
candidate. 

• BBCp : All other BB events with a misrecon- 
structed D candidate. 

• BBCd '■ Other BB events with a correctly recon- 
structed D — > ir + ir~ir° decay. 

• qqp : Continuum e + e _ — > qq events with a misre- 
constructed D candidate. 

• qqo : Continuum e + e~ — > qq events with a cor- 
rectly reconstructed D — > ir + ir~ir° decay. 

The Cabibbo-favored decay chain B~ —> D°ir~, D° —> 
K~ir + ir°, which has the same final state particles as 
our signal decay, does not contribute significantly to the 
background, since it is suppressed by the particle identi- 
fication and mo cuts 

The majority of background events are of the qqp type. 
The combination of DttX, DKX, and BBCp events con- 
stitutes the second largest background. In order to sup- 
press these backgrounds, we have developed two neural 
networks, each of which combines several input variables 
that provide separation between signal and background. 
The first neural network variable q is computed from in- 
put variables that provide separation between continuum 
and BB events. The second variable d combines input 
variables that separate correctly reconstructed tt° and D° 
candidates from misreconstructed ones. It provides sep- 
aration between signal and all misreconstructed- D back- 
ground. 

The input variables for q are (1) the cosine of the CM 
angle between pe and the beams; (2) | cos#r|; (3-4) the 
zeroth and second Legendre moments of the momentum 
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flow of the ROE about the CM thrust axis of the B can- 
didate daughters; (5) log of the distance along the beam 
direction between the reconstructed B vertex and the 
vertex of the ROE, computed as in Ref. EH; (6) log of 
the distance of closest approach between the kaon track 
and the D decay vertex, which is calculated from the tt + 
and 7T _ track parameters; (7) an integer variable calcu- 
lated from the probability that the ROE contains a B°, 
determined using the lepton flavor tagging algorithm of 
Ref. |H|. 

The input variables for d are (1) the invariant mass 
of the 7T° candidate; (2) the tt° momentum in the lab 
frame; (3) cosine of the tt° decay angle 9^o , defined as 
the angle between the momentum and the momentum 
of the CM frame, calculated in the ir° rest frame; (4) the 
invariant mass m/j of 7r°, where 7r° is the ir° candidate 
reconstructed from the jh and any additional photon in 
the event except 7 S , chosen such that m/, is closest to the 
nominal ir° mass; (5) m s , calculated analogously to m^, 
but with 7 S instead of jh', (6-7) cosines of the decay angles 
of the 7r° and the 7r°, calculated analogously to w o; (8) 
cosine of the angle between pe and the thrust axis of 
the 7r + 7r _ 7r° final state, calculated in the tt + tt~tt° rest 
frame; (9) cosine of the angle between the D candidate 
momentum and the line connecting the B and D decay 
vertices. 

The q and d distributions of simulated signal and back- 
ground events are shown in Fig.^ All events are required 
to satisfy the conditions q > 0.1, d > 0.1, in order to re- 
duce the background and suppress correlations between 
the variables used in the fit described below. The final 
signal reconstruction efficiency is 10.5%. 




FIG. 1: (a) Distribution of the neural network variable q for 
continuum (hatched) and signal simulated events. The BB 
background distribution is similar to that of signal, (b) Dis- 
tribution of the neural network variable d for BB background 
(hatched) and signal simulated events. The continuum back- 
ground distribution is similar to that of the BB background. 
All histograms are normalized to unit area. 



We perform a maximum likelihood fit to measure the 
number and the decay rate asymmetry of signal events 
in the on-resonance data sample, using the variables AE, 
q, and d. The variable rriESi which is commonly used as 



a fit variable in B decay analyses, is not included in the 
fit. Studies with simulated events indicate that correla- 
tions of m es with other fit variables in the distributions 
of DttX, DKX, and BBCp background events lead to 
a bias in the measured signal yield, unless the correla- 
tions are modeled correctly. Such modeling complicates 
the analysis procedure, increases the dependence on the 
simulation, and incurs additional systematic errors. By 
excluding rriES from the fit, We give up some statisti- 
cal precision in order to make the analysis more robust. 
Correlations between the AE, q, and d distributions for 
the different event types are at the few percent level in 
the worst cases, and ignoring them in fits to simulated 
events does not result in significant biases. 

The probability density function (PDF) for the fit is 

V = ~Y j N{P t {AE,q,d), (1) 

where the subscript t corresponds to one of the ten event 
types listed above, TV* is the number of events of type t, 
and r\ = J2t Nt- The PDF Vt for events of type t is a 
product of the form 

V t (AE,q,d) = £ t (AE)Q t (q)V t (d). (2) 

The functions £ BBCp (AE), £ qqD (AE), and £ qqp (AE) 
are parameterized as second order polynomials, and all 
other £ t (AE) functions are the sum of a Gaussian and a 
second order polynomial. The parameters of these func- 
tions are obtained from fits to simulated events. The 
PDFs Qt(q) and V t (d) are 15-bin histograms obtained 
from simulated events. 

To extract the signal yield and asymmetry, we mini- 
mize the log of the extended likelihood 

^Hrll^. (3) 

i=l 

Six parameters are floating in the fit. These are the 
event yields N DKd , N D7TD , N qqp , and N BBp = N DKX + 
NdttX + NbbC p , the ratio R D7T x = N d ^ x /N B b p , 
and the decay rate asymmetry A = (N BKd — 
n dk d )I( n dk d + n dk d )i where the superscript indi- 
cates the charge of the kaon. Five ratios of event yields 
are obtained from the simulation and are not varied 
in the fit. From these ratios we obtain the five pa- 
rameters N DK x = 0.21 Ndkx, N DlTp = 0.171 N D7TD , 
N BBCd = 0.0089 N BBlB , N qqD = 0.0136 N qqp , and 
Ndk p — 0.1614 N B k d ■ All fixed parameters are later 
varied to evaluate systematic errors, as described below. 

The results of the fit are summarized in Table [I] We 
observe N B k d = 133 ± 23 signal events and the decay 
rate asymmetry A — 0.02 ±0.16, where the errors are 
statistical only. The corresponding branching fraction 
is B(B~ -> D n+W -„ K~) = (5.5 ± 1.0) x lO" 6 . The 
statistical significance of the signal observation, obtained 
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from a scan of the likelihood as a function of the signal 
yield, is 5.7 standard deviations. 

The fit parameter most correlated with the sig- 
nal yield is Ng Bp , with correlation matrix element 
p{Ndk d ,Nbbp) = —0.33. The largest correlation ma- 
trix element for the asymmetry is p(A, Ndk d ) — —0.036. 
Projections of the data and the fit function onto the fit 
variables are shown in Fig. for events with a high like- 
lihood of being signal and for the entire data sample. 



TABLE I: Results of the data fit. Errors are statistical only. 



Parameter 


Value 


n DKd 

A 


133 ± 23 
0.02 ± 0.16 


Nqqp 

N B Bp 

RdttX 


43 ± 16 
1454 ± 53 
806 ± 54 
0.82 ± 0.11 



The systematic uncertainties in the signal branching 
fraction and asymmetry measurements are summarized 
in Table ITT1 We describe briefly the procedures used for 
their evaluation: (1) The statistical errors in the simu- 
lated samples used to obtain the shapes of £ t ( AE) , Qt(q), 
and T> t (d) are propagated to the final fit. (2) The value 
of Ndkx is varied by ±25%, determined from the uncer- 
tainties on the decay modes contributing to the DKX 
background Q. The parameters N DlTp , N B BCm N qqD1 
and Ndk p are varied by ±50%, which is estimated to 
be very conservative, given the level of data-simulation 
agreement. (3) We evaluate the effect of possible dif- 
ferences between the event distributions in the data and 
the simulation by studying events in the tues sideband 
5.23 < rriEs < 5.26 GeV/c 2 , as well as events produced in 
the copious decay mode B~ — > D°tt~, D° — > K~tt + tt°. 
(4) The AE and d distributions of DKo events are 
slightly correlated, and this correlation is ignored in the 
PDF. To evaluate the uncertainty due to this, we re- 
peat the fit with T>oK D (d) taken from simulated DKjj 
events in different AE bins. (5) We consider the effect of 
a possible contribution of charmless B~ — > K~ tt + it~ 7r° 
events, assuming a branching fraction of 6 x 10~ 5 . (6) 
The uncertainty in the contribution of a non-resonant 
component to the D — > 7r + 7r~7r° decay || is propagated 
to the signal efficiency. (7) We account for the possi- 
bility of charge-dependence in the track reconstruction 
efficiency. (8-9) We assign a reconstruction efficiency un- 
certainty of 1.4% per charged track and 3.5% for the 7r°. 
(10) We account for the uncertainty in the number of BB 
events produced by PEP-II and (11) the uncertainty in 
the efficiency of the particle identification requirements 
applied to the data sample. 

Additional cross-checks are performed to verify the va- 
lidity of our results. We compare the fit variable distri- 
butions of the data with those of simulated events in the 




in 

§140 

2-120 

"100 
c 

0) 

> 

HI 



-0.06 -0.04 -0.02 0.02 0.04 0.06 
A E(GeV) 





BaBar 















0.06 -0.04 -0.02 0.02 0.04 0.06 
A E(GeV) 




J.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

d 



1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

d 



FIG. 2: Projections of the data (data points) and fit func- 
tion onto the (a) AE, (b) q, and (c) d axes. The curves or 
histograms in each plot show, from bottom to top, the cumu- 
lative contributions of the continuum, BB background, and 
signal components of the fit function. For each of the vari- 
ables plotted, a tight event selection is applied on the other 
two variables with a signal efficiency of about 50%, in order 
to increase the signal-to-background ratio in these plots. Fig- 
ures (d), (e), and (f) show the same projections for the entire 
data sample. 



TABLE II: Fractional systematic error in the signal branching 
fraction B and absolute error in the asymmetry A. 



No. Source 


Error (%) 




(J B /B 


OA 


(1) Simulated sample statistics 


7.9 


1.8 


(2) Variation of fixed yields 


6.2 


0.25 


(3) Data-simulation shape comparison 


5.8 


1.6 


(4) AE — d correlations in Vdk d 


1.9 


0.39 


(5) Charmless branching fraction 


0.85 


0.093 


(6) Dalitz plot distribution 


0.33 




(7) Detector asymmetry 




0.90 


(8) Track reconstruction efficiency 


4.2 




(9) 7T° efficiency 


3.5 




(10) Number of BB events produced 


1.1 




(11) Particle ID efficiency 


1.0 




Total 


13 


2.6 
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AE sideband 90 < AE < 140 MeV. The fit variable dis- 
tributions of simulated continuum events are validated 
against the off-resonance data. The simulated distribu- 
tions of B~ — ► D°tt~ , D° — » 7r + 7T _ 7r° events are com- 
pared with their distributions in the data. We verify 
the signal efficiency by measuring the branching fraction 
B(B~ -> D°tt-) using D° decays to D° -> if-vr+vr 
and Z?° — > 7r + 7r~7r°. The simulated distributions of the 
q and d input variables are compared with the distri- 
butions in the data. In all cases, good agreement be- 
tween simulation and data is observed. No significant 
excess of signal events is found in a fit to data events 
in the m D sidebands 1.775 < m D < 1.800 GeV/c 2 
and 1.920 < m D < 1.955 GeV/c 2 . We conduct fits 
to event samples containing simulated signal and back- 
ground events and find no significant biases in all the fit 
variables. Fits to parameterized experiments generated 
with the parameter values obtained in the data fit are 
unbiased, and their distributions of fit parameter errors 
and maximum likelihood are consistent with those of the 
data fit. 

In summary, using a sample of 229 ± 2.5 million 
e + e~ — > BB events we observe 133 ± 23 events in the 
decay chain B~ — > D^+^-^oK - , where the 7r + 7r~7r° fi- 
nal state excludes the CP-eigenstate K®tt . We extract 
the branching fraction and decay rate asymmetry 



B(B- 



-) = (5.5 ± 1.0 ±0.7) x 10~ 6 , 
A = 0.02 ±0.16 ±0.03, (4) 



where the first errors are statistical and the second are 
systematic. The level of background suppression we 
achieve is critical for using this mode to measure 7. The 
remaining background doubles the statistical error on 7 
with respect to the no-background case. 
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